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Mathematical  models  of  an  indirect  internal  reforming  solid  oxide  fuel  cell  (IIR-SOFC)  fed  by  four  dif¬ 
ferent  primary  fuels,  i.e.,  methane,  biogas,  methanol  and  ethanol,  are  developed  based  on  steady-state, 
heterogeneous,  two-dimensional  and  tubular-design  SOFC  models.  The  effect  of  fuel  type  on  the  thermal 
coupling  between  internal  endothermic  reforming  with  exothermic  electrochemical  reactions  and  sys¬ 
tem  performance  are  determined.  The  simulation  reveals  that  an  IIR-SOFC  fuelled  by  methanol  provides 
the  smoothest  temperature  gradient  with  high  electrochemical  efficiency.  Furthermore,  the  content  of 
CO2  in  biogas  plays  an  important  role  on  system  performance  since  electrical  efficiency  is  improved  by 
the  removal  of  some  CO2  from  biogas  but  a  larger  temperature  gradient  is  expected. 

Sensitivity  analysis  of  three  parameters,  namely,  a  operating  pressure,  inlet  steam  to  carbon  (S:C)  ratio 
and  flow  direction  is  then  performed.  By  increasing  the  operating  pressure  up  to  lObar,  the  system  effi¬ 
ciency  increases  and  the  temperature  gradient  can  be  minimized.  The  use  of  a  high  inlet  S:C  ratio  reduces 
the  cooling  spot  at  the  entrance  of  reformer  channel  but  the  electrical  efficiency  is  considerably  decreased. 
An  IIR-SOFC  with  a  counter-flow  pattern  (as  based  case)  is  compared  with  that  with  co-flow  pattern  (co¬ 
flow  of  air  and  fuel  streams  through  fuel  cell).  The  IIR-SOFC  with  co-flow  pattern  provides  higher  voltage 
and  a  smoother  temperature  gradient  along  the  system  due  to  superior  matching  between  heat  supplied 
from  electrochemical  reaction  and  heat  required  for  steam  reforming  reaction;  thus  it  is  expected  to  be 
a  better  option  for  practical  applications. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

During  the  past  two  decades,  several  types  of  fuel  cell  have 
been  continuously  developed  and  modified.  Among  them,  the  solid 
oxide  fuel  cell  (SOFC)  has  attracted  much  interest  for  station¬ 
ary  power  generation  since  it  provides  high  efficiency  with  low 
pollutant  emission  [1  ].  Generally,  SOFC  is  operated  at  high  temper¬ 
atures  (1073-1473  K),  thus  hydrocarbon  fuels  (e.g.,  methane)  can 
be  applied  as  a  primary  fuel  when  the  system  is  operated  with 
internal  reforming  (IR-SOFC).  This  operation  uses  the  heat  gener¬ 
ated  from  the  exothermic  electrochemical  reaction  to  conduct  the 
endothermic  (steam)  reforming  of  hydrocarbons  [2-5]. 

Theoretically,  there  are  two  main  designs  of  a  IR-SOFC 
namely,  direct  internal  reforming  (DIR-SOFC)  and  indirect  inter¬ 
nal  reforming  (IIR-SOFC).  In  the  DIR-SOFC  approach,  together  with 
electrochemical  reaction,  the  reforming  reaction  occurs  simultane¬ 
ously  at  the  anode  side  of  SOFC.  Thus,  a  high  heat  transfer  rate  can 
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be  achieved  from  this  operation;  nevertheless,  the  anode  material 
must  be  optimized  for  both  reactions  and  can  be  easily  poisoned 
by  deposition  of  carbon  produced  from  the  reforming  of  hydro¬ 
carbons.  In  an  IIR-SOFC,  the  endothermic  reforming  reaction  takes 
place  at  a  reformer,  which  is  in  close  thermal  contact  with  the  anode 
side  of  the  SOFC  where  exothermic  electrochemical  reaction  occurs. 
The  great  advantage  of  an  IIR-SOFC  is  that  the  released  heat  can  be 
utilized  efficiently  and  the  anode  material  is  not  subject  to  coke 
deposition  occurring  from  the  cracking  and  reforming  of  hydro¬ 
carbons,  thus  the  primary  fuels  for  IIR-SOFC  are  more  flexible  than 
DIR-SOFC.  It  is  noted  that  the  main  drawback  of  IIR-SOFC  is  the  mis¬ 
match  between  the  rates  of  endothermic  and  exothermic  reactions. 
This  problem  leads  to  significant  local  temperature  reduction,  par¬ 
ticularly  near  the  entrance  of  the  internal  reformer  [3,4,6],  which 
can  result  in  mechanical  failure  due  to  thermally  induced  stresses. 
Importantly,  the  temperature  gradient  along  the  SOFC  system  is 
affected  by  several  operating  parameters,  namely,  inlet  tempera¬ 
ture,  pressure  and  drawn  current  [4,6],  as  well  as  the  type  of  primary 
inlet  fuel. 

Focusing  on  the  primary  fuel  selection,  methane  (or  natural  gas) 
is  accepted  to  be  the  most  applicable  fuel  for  IR-SOFC  due  to  its 
abundance  and  well-developed  supply  infrastructure.  Among  other 
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Nomenclature 

Aact  area(m2) 

Aact  active  area  =  7l{dr~2 rcat)L  (m: 2 ) 

7l{d2r-{dr-2rcat)2l) 

CP  specific  heat  of  the  gas  streams  (kj  mol-1  K-1 ) 
q  concentration  of  component  i  (mol  m-3 ) 

Dij  binary  diffusion  of  component  i  and  j  (m2  s-1 ) 

Deimix  effective  molecular  diffusivity  (m2  s-1 ) 

d  diameter  (m) 

Eff  efficiency 

F  Faraday’s  constant  (96,484  C  mol-1 ) 

Fr,s  gray-body  transfer  factor  from  reformer  surface  to 
solid  oxide  surface 
AH  heat  of  reaction  (kj  mol-1 ) 

Jo  Exchange  current  density  (A  m-2 ) 

J  current  density  (A  m-2 ) 

kp  permeability 

LHVi  lower  heating  value  of  component  i 

Mi  molecular  weight  of  component  i  (kg  kj-1 ) 

n/  number  of  mole  of  component  i 

p°  standard  partial  pressure  (Pa) 

Pi  partial  pressure  of  species  i 

Psofc  power  density  ( W  m-2 ) 

Qe/ec  local  generated  heat  ( W  nrr 2 ) 

Rk  rate  of  reaction  k 

Rg  universal  gas  constant:  8.414  kj  mol-1  K-1 . 

R0hm  ohmic  resistant  ( k£2  nrr 2 ) 

T  temperature  (K) 

v  fluid  velocity  (ms-1) 

yi  the  mole  fraction  of  gas 

Greek  letters 

p  density  (kg  m-3) 

sm  emittance 

fi  effect  of  viscosity  (Pa  s) 

y  special  Fuller  et  al.  diffusion  volume  [  13  ] 

ofa.c  charge  transfer  coefficient  of  anode  and  cathode 
a  Stefan-Boltzmann  coefficient 

s  porosity 

r  tortuosity 

A  thermal  conductivity  (kj  m-1  s-1 1<-1 ) 

r\  voltage  drop  (V) 

r  thickness  (m) 

Subscripts 

a  air  channel 

i  component  (methane,  water,  hydrogen,  etc.) 

j  reaction  (SRM,  WGS,  etc.) 

/  fuel  channel 

p  particle 

r  reformer 

s  solid  oxide  fuel  cell 

Act  activation  losses 

Cat  catalyst 

Cell  cell  stack 

Con  concentration  losses 

Elec  electrical 

Thermal  thermal 


potential  primary  fuels,  ethanol  and  methanol  are  also  interest¬ 
ing  candidates  due  to  their  ready  availability,  high  specific  energy, 
sulfur-free  content  and  storage  transportation  convenience  [5]; 
moreover,  they  can  be  produced  renewably  from  both  chemical  and 


biological  processes.  In  addition  to  these  three  types  of  fuel,  biogas 
is  another  attractive  candidate— particularly  for  agricultural  coun¬ 
tries  since  it  is  a  promising  renewable  fuel  that  can  be  produced 
biologically  from  the  waste  residues.  Normally,  after  desulfuriza¬ 
tion,  biogas  is  composed  of  approximately  60-65  vol.%  of  methane 
and  35-40  vol.%  ofC02. 

In  the  present  study,  a  set  of  mathematical  models  has  been 
developed  to  predict  and  compare  the  behaviour  of  a  tubular 
IIR-SOFC  fuelled  by  methane,  methanol,  ethanol,  and  biogas.  The 
models  are  coded  in  the  COMSOL®  program  within  2D  axisymmet- 
ric  application.  From  the  simulation,  the  performance  of  an  IIR-SOFC 
(i.e.,  electrical  efficiency  and  temperature  gradient  along  the  sys¬ 
tem)  fuelled  by  these  four  types  of  primary  fuels  are  determined; 
in  addition,  the  effects  of  inlet  steam  to  carbon  (S:C)  ratio,  operating 
pressure  and  gas  flow  patterns  on  the  performance  of  the  IIR-SOFC 
are  also  investigated.  Finally,  a  suitable  type  of  primary  fuel  and 
suitable  operating  conditions  for  an  IIR-SOFC  system  are  identified. 

2.  Mathematical  models 

2.1.  Kinetic  equations  for  steam  reforming  of  primary  fuels 

A  steam  reforming  reaction  was  chosen  to  convert  methane, 
methanol,  ethanol  and  biogas  to  hydrogen-rich  gas.  The  kinetic 
equations  for  each  reaction  were  based  on  previous  reports  in  the 
literature  [7-10].  For  methane  reforming,  the  well-known  kinetic 
models  and  reaction  mechanisms  proposed  by  Xu  and  Froment 
[7]  were  applied.  For  methanol  reforming,  the  rate  expressions 
reported  by  Peppley  et  al.  [8]  were  selected,  whereas  the  simulation 
of  ethanol  reforming  was  based  on  the  rate  expressions  proposed  by 
Sahoo  et  al.  [10].  For  biogas,  the  rate  expressions  of  methane  steam 
reforming  with  coupling  of  carbon  dioxide  reforming  reported  by 
Abashar  [9]  were  applied. 

2.2.  Model  geometry  and  assumptions 

A  tubular-designed  IIR-SOFC  was  selected,  as  shown  schemati¬ 
cally  in  Fig.  1.  According  to  this  configuration,  primary  fuels  (e.g., 
methane,  methanol,  ethanol,  or  biogas)  and  steam  are  converted 
to  hydrogen-rich  gas  in  an  internal  catalytic  packed-bed  reformer 
before  being  introduced  to  the  fuel  channel  of  the  SOFC.  Simul¬ 
taneously,  air  is  fed  in  an  opposite  flow  direction  through  the  air 
channel.  It  is  noted  that  all  dimensions  and  physical  properties  of 
SOFC  system  in  the  present  work,  which  are  summarized  in  Table  1 , 
are  based  on  the  work  of  Aguiar  et  al.  [4]. 

The  mathematical  model  was  developed  as  the  smallest,  single, 
unit  cell  taking  into  account  the  effect  of  temperature  on  gas  dis¬ 
tribution,  reactant  conversion  and  charge  transfer.  Two  main  sets 
of  equations,  namely,  for  the  internal  reformer  and  the  SOFC  stack 
(including  electrodes  and  solid  electrolyte)  were  applied  to  pre¬ 
dict  the  concentration  and  temperature  gradients  along  the  system 
based  on  the  following  assumptions: 


M -  Dimensionless  Axial  - ► 


Fig.  1.  Schematic  diagram  of  tubular-designed  SOFC  with  indirect  internal  reforming 
operation. 
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Table  1 

Dimensional  parameter  for  modelling. 


Parameters 

Value 

Reformer  diameter 

0.6  cm 

Gap  between  reformer  to  fuel  channel 

0.24  cm 

Anode  thickness 

0.07  cm 

Electrolyte  thickness 

20  pan 

Cathode  thickness 

50  pan 

Gap  between  cathode  to  air  channel 

0.3  cm 

Length 

60  cm 

Anode 

Porosity  (£) 

0.35 

Tortuosity  (r) 

4.80 

Pore  radius  (rp) 

0.20  [Jim 

Particle  diameter  (dp) 

1.00  [Jim 

Cathode 

Porosity  (e) 

0.35 

Tortuosity  (r) 

4.00 

Pore  radius  (rp) 

0.25  [Jim 

Particle  diameter  (dp) 

1.25  [Jim 

•  Each  section  is  considered  to  be  under  non-isothermal,  steady- 
state  conditions. 

•  Ideal  gas  behaviour  is  applied  for  all  gas  components. 

•  Diffusion  in  the  gas  phase  and  pressure  drop  in  the  SOFC  stack 
and  reformer  are  ignored. 

•  The  operating  cell  potential  is  constant  along  the  cell  coordinate. 

•  Final  fuel  utilization  is  fixed  constantly  at  80%. 

By  adopting  a  non-isothermal  state,  the  calculated  gas  proper¬ 
ties  are  derived  as  the  function  of  temperature;  thus,  momentum, 
mass  and  energy  balances  are  integrated  under  the  conditions  of 
interest.  According  to  the  tubular  IIR-SOFC  design  (as  shown  in 
Fig.  1 ),  the  system  structure  consists  of  two  main  parts,  i.e.,  porous 
and  bulk  gas  channels.  The  transport  behaviour  of  both  parts  is 
treated  differently,  as  explained  in  the  following  sections. 


2.2.1.  Gas  transport 

Gas  transport  behaviour  in  the  porous  media,  reformer  and  cell 
layer  was  calculated  by  means  of  the  Brinkmen  equation  (Eq.  (1)), 
whereas  the  Navier-Stokes  equation  (Eq.  (2))  was  applied  to  predict 
the  momentum  transfer  in  gas  channels  [11  ]. 

Vp  =  -j^v  +  £LV2v  (1) 

Kp 

v  V(pv)  =  -Vp  +  /xV2v  (2) 


where  v  is  the  fluid  velocity;  p  is  the  density  (kg  m-3 );  p  is  the  pres¬ 
sure  (Pa);  /x  is  the  effect  of  viscosity  (as  function  of  temperature), 
and  kp  is  the  permeability.  In  this  system,  the  diffusion  properties  of 
several  gas  species  were  investigated  by  using  the  binary  diffusion 
flux,  as  given  by  [12]: 


DiJ  = 


(0. 00143  )T175 


(3) 


where  D/j  is  the  binary  diffusion  of  component  i  (molm-1  s-1);  T 
is  the  temperature  (I<);  p  is  the  pressure  (Pa);  M2J-  =  2/(1 /M,  +  1/M,); 
Mi  is  the  molecular  weight  of  component  i;  y  is  a  special  diffusion 
volume  (reported  by  Fuller  et  al.  [13]);  p2  is  the  partial  pressure  of 
component  i.  In  the  case  of  transportation  through  porous  media, 
the  diffusion  behaviour  is  corrected  by  taking  into  account  the 
porosity,  s,  and  the  tortuosity,  r,  namely,  the  effective  diffusivity 
coefficiency  (m2  s-1 ),  Deimix: 


De  . 

i,mix 


(4) 


Generally,  gas  diffusion  through  electrode  layers  can  be 
explained  by  the  coupling  of  two  diffusion  mechanisms,  namely, 


molecular  and  Knudsen  diffusions.  This  is  known  as  the  dusty  gas 
model  (DGM)  [14].  If  the  mean  free  path  of  gas  molecular  species 
is  larger  than  the  pore  diameter,  Knudsen  diffusion  is  the  main 
mechanism;  otherwise,  molecular  diffusion  will  dominate. 


2.2.2.  Ionic  transport 

The  electricity  generation  was  calculated  from  the  relation 
between  polarization  and  the  activated  area  of  SOFC.  Activation, 
concentration  and  ohmic  polarization  are  known  as  major  losses  for 
high-temperature  fuel  cells  like  SOFCs.  Activation  loss  arises  from 
an  activation  barrier  to  the  electrochemical  reaction  at  the  elec¬ 
trode.  In  this  work,  the  activation  polarization,  pact,  was  computed 
by  [15]: 

^  -  exp 

where  J2  is  the  current  density  (Am-2);  J0  is  the  exchange  current 
density  (Am-2);  Rg  is  the  universal  gas  constant  (kjmol-1 1<-1);  F 
is  the  Faraday  constant.  All  relevant  parameters  were  reported  by 
Zhu  et  al.  [16].  The  concentration  polarization  occurs  due  to  the 
resistance  to  gas  diffusion  through  the  porous  media.  Generally, 
gas  diffusion  behaviour  can  be  predicted  by  three  mathematical 
models:  Fields  model,  Dusty  gas  model  (DGM),  and  Stefan-Maxwell 
model;  the  DGM  was  chosen  in  the  present  work.  Ohmic  polar¬ 
ization  arises  from  the  ion  transport  across  the  cell,  which  mainly 
depends  on  the  ionic  conductivity  of  the  SOFC  material.  By  applying 
ohm’s  law,  the  relation  of  ohmic  polarization  and  material  resistiv¬ 
ity  can  be  determined.  It  is  noted  that,  in  the  present  work,  the 
material  properties  reported  by  Zhu  et  al.  [16]  were  used.  This  volt¬ 
age  drop  was  computed  by  the  simplified  equation: 


-(1  -  ac)Ftiact'j 


RgT 


(5) 


Ji  —  Jo 


exp 


(  OiaFpact 


*1  ohmic  —J-Rohm  (6) 

where  p0hmic  is  the  voltage  drop  caused  by  ohmic  losses  (V);  J  is  the 
current  density  (Am-2),  and  Rohm  is  the  ohmic  resistant  per  unit 
area  (k£2  m2).  It  is  noted  that  although  both  H2  and  CO  can  be  elec- 
trochemically  consumed  in  a  SOFC,  it  has  been  reported  that  the 


Fig.  2.  I-V  curve  validation  between  simulation  results  of  present  work  (blank  sym¬ 
bols)  with  those  from  Lin  et  al.  [16]  (solid  symbols). 


Table  2 

Thermal  properties  of  each  fuel  under  isothermal  condition  (1173  K  and  1  bar). 


Fuel 

Thermal  properties 

LHV 

Thermal  efficiency  (%) 

Methane 

744.7 

49.4 

Biogas 

466.8 

39.5 

Methanol 

631.5 

76.3 

Ethanol 

1194.3 

69.7 
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rate  of  electricity  generation  from  CO  is  three  times  lower  than  that 
of  H2  [17];  thus  the  rate  of  CO  oxidation  in  SOFC  has  been  neglected 
in  this  study.  At  the  anode  [electrolyte  interface,  H2  is  electrochem- 
ically  oxidized  whereas  02  is  consumed  at  the  cathode  [electrolyte 
interface.  The  conversion  rate  of  H2,  Re/ec?H2,  and  02,  Reiec,o2,  at  the 
indicated  boundary  was  simulated  from  Eqs.  (7)  and  (8),  respec¬ 
tively: 


Kelec,H2  ~  2f 


(7) 


P 

Kelec,02  ~  4j7 


(8) 


The  reliability  of  these  electrochemical  equations,  coded  in 
COMSOL®  program,  was  validated  by  comparison  with  the  results 
from  the  work  of  Lin  et  al.  [18],  as  shown  in  Fig.  2.  Clearly,  the  results 
are  in  good  agreement. 


2.2.3.  Heat  transfer 

Conduction  heat  transfer  along  stack  materials,  convection  heat 
transfer  from  fluid  flow  though  the  system  and  radiation  between 
the  reformer  and  the  solid  electrolyte  were  considered  to  predict 
the  temperature  distribution  of  the  IIR-SOFC  system.  In  the  bulk 
gas  channel,  both  conduction  and  convection  heat  transfer  were 
calculated  by: 

V  •  ( vpcpT )  -  V  •  (A.  •  VT)  =  0  (9) 

where  v  is  the  fluid  velocity  (ms-1);  p  is  the  gas  density  (kgm-3); 
cp  is  the  specific  heat  (kjkg-1  K-1);  A  is  the  thermal  conductiv¬ 
ity  (kjm-1  s-1  K-1);  T  is  the  temperature  (I<).  It  is  noted  that  the 


influence  of  the  radiation  between  the  reformer  and  the  solid  elec¬ 
trolyte  was  also  taken  into  account  by  considering  the  heat  required 
for  steam  reforming  in  the  reformer  (Eq.  (10))  and  the  heat  of 
electrochemical  reaction  at  the  anode  [electrolyte  layer  (Eq.  (11)) 
[19]. 


V  •  ( VpCpT )  V  •  (A  •  VT)  +  (  AH refomi  )Rreform 


_ gCrr4-7?) 

(1  —  £m / 6mA)r  +  l/ArFrjS  +  (1 


6m/ 6mA)s 


(10) 


V  •  ( vpcpT )  -  V  •  (A  •  VT)  +  Qe/eC 


_ v(T?-Tr4) _ 

(1  —  £m/ £mA)s  +  1  /  AsFr,s  +  (1  —  6m/ 6mA)r 


=  0 


(11) 


where  Qeiec  =  (-AHeiec)ReieC  -  Psofc  is  the  local  generated  heat 
(Wm-2);  (-A Hk)  is  the  heat  of  reaction  k;  A  is  the  area  (m2);  Fr,s 
is  the  gray-body  transfer  factor  from  the  reformer  surface  to  the 
solid  oxide  surface;  a  is  the  Stefan-Boltzmann  coefficient;  £m  is  the 
emittance;  Psofc  is  the  local  power  density  of  the  SOFC  (W  m-2). 


2.3.  Calculation  of  thermal  and  electrical  efficiencies 

Base  on  the  assumption  that  most  of  the  generated  hydrogen 
was  combusted  to  supply  heat  for  the  steam  reforming  reaction, 
some  heat  generated  was  used  as  the  heat  source.  The  thermal  effi¬ 
ciency,  Fff thermal ,  of  an  IIR-SOFC  fuelled  by  different  primary  fuels 
was  defined  as  the  ratio  between  the  heating  value  of  hydrogen 


Fig.  3.  Concentration  profiles  in  reformer  (black  symbols)  and  fuel  channel  (white  symbols)  of  IIR-SOFC  fuelled  by  (a)  methane,  (b)  biogas,  (c)  methanol  and  (d)  ethanol  (at 
base  case  conditions). 
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leftover  from  combustion  with  that  of  the  inlet  fuel,  i.e.: 

nft'LHV*.  1173  K  -  E  (-^reform) 


Eff thermal  ~ 


En?LHVu 


(12) 


i,1173  I< 


where  nz  is  number  of  moles  for  component  i,  and  LHVi  n 73K  is  the 
lower  heating  value  of  component  i  at  1173  K.  In  order  to  investigate 
the  electrical  efficiency,  Effdec,  outlet  reformats  were  fed  into  SOFC 
and  the  electrical  efficiency  was  defined  as: 


Eff dec  ~ 


PsOFC-Aact 


(13) 


EyTEHvun73K 

where  Aact  is  the  active  area  (m2)  and  yjfje/  is  the  mole  fraction  of 
the  primary  fuel. 


3.  Result  and  discussion 

3. 1.  Thermodynamic  property  comparison 

Before  undertaking  a  thermodynamic  analysis  of  an  IIR-SOFC 
fuelled  by  different  primary  fuels,  the  thermodynamic  properties 
of  each  fuel,  i.e.,  methane,  methanol,  ethanol,  and  biogas  (60%  CFI4 
and  40%  CO2 )  were  identified  in  terms  of  lower  heating  value  (LFIV) 
and  thermal  efficiency.  According  to  the  thermal  efficiency  calcu¬ 
lation,  a  steady-state  operation  was  assumed  and  the  inlet  fuel 
was  kinetically  converted  to  hydrogen  in  the  packed-bed  reformer, 
which  was  eventually  combusted  to  predict  the  thermal  efficiency. 
It  is  noted  that  the  simulation  was  carried  out  at  1173  K  and  1  bar 
with  the  inlet  steam  to  carbon  (S/C)  ratio  of  2.0.  Importantly,  the 
limitation  of  gases  transport  through  the  packed  catalyst  was  also 
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Fig.  4.  Results  from  COMSOL®  represented  temperature  gradients  in  reformer  and  fuel  channel  of  IIR-SOFC  fuelled  by  (a)  methane,  (b)  biogas,  (c)  methanol  and  (d)  ethanol 
(at  base  case  conditions). 
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Fig.  5.  Temperature  profiles  of  reformer,  fuel  and  air  channels  (Tr,  Tf  and  Ta,  respectively)  for  IIR-SOFC  fuelled  by  (a)  methane,  (b)  biogas,  (c)  methanol  and  (d)  ethanol  (at 
base  case  conditions). 


taken  into  account.  The  dimensions  of  the  system  configuration  are 
listed  in  Table  1,  while  the  LHV  and  thermal  efficiencies  of  each 
primary  fuel  are  compared  in  Table  2.  Among  these  primary  fuels, 
alcohol-based  fuels  (methanol  and  ethanol)  provide  significantly 
higher  thermal  efficiency  than  methane-based  fuels.  Furthermore, 
the  LHV  of  ethanol  is  considerably  higher  than  those  of  other  fuels, 
while  that  of  methanol  is  close  to  that  of  methane.  Hence,  alcohol 
would  be  a  promising  candidate  fuel  for  IIR-SOFC,  IIR-SOFCs.  The 
performance  of  an  IIR-SOFC  fuelled  by  these  primary  fuels  in  terms 
of  electrical  efficiency  and  temperature  gradient  along  the  system 
was  then  compared  as  presented  in  the  following  sections. 

3.2.  Modelling  of  IIR-SOFC  fuelled  by  different  fuels  as  based  case 

The  non-isothermal  model  of  IIR-SOFC  was  first  simulated  at 
1  bar  by  introducing  inlet  fuels  and  steam  at  a  S:C  ratio  of  2.0.  A 
SOFC  load  voltage  of  0.7  V  and  a  fuel  utilization  of  80%  were  applied. 
Under  these  conditions,  characteristic  results,  i.e.,  primary  fuel  con¬ 
version,  product  gas  distribution  and  temperature  gradient  along 
the  internal  reformer  and  SOFC  channels  (both  fuel  and  air  chan¬ 
nels)  are  shown  in  Figs.  3-5.  The  mole  fraction  profiles  of  primary 
fuels,  hydrogen,  steam,  carbon  monoxide,  and  carbon  dioxide  in 
reformer  and  fuel  channel  of  SOFC,  are  presented  in  Fig.  3(a)-(d). 
Among  these  primary  fuels,  ethanol  provides  the  slowest  conver¬ 
sion  rate  along  the  reformer  channel.  At  the  outlet  of  the  internal 
reformer,  the  product  gas  flows  backward  into  fuel  channel  where 
hydrogen  and  carbon  monoxide  are  electrochemically  converted  to 
steam  and  carbon  dioxide. 

Fig.  4(a)-(d)  shows  the  results  from  COMSOL  software  repre¬ 
senting  temperature  distributions  in  reformer,  fuel  and  air  channels 


for  the  IIR-SOFC  system  fuelled  by  different  fuels,  while  Fig.  5(a)-(d) 
summarizes  the  temperature  profiles  in  each  channel.  It  can  be 
seen  that  the  heat  generated  from  fuel  cell  side  is  sufficient  for  the 
internal  reforming  operation;  nevertheless,  a  noticeable  tempera¬ 
ture  gradient  associated  with  indirect  internal  reforming  operation 
is  also  observed  for  all  types  of  fuel  (as  also  reported  by  other 
researchers  [4,6] ).  Among  them,  IIR-SOFC  fuelled  by  methane  shows 
the  largest  cooling  spot  at  the  first  half  of  the  reformer  channel 
(reformer  temperature  decreases  from  1173  to  1000  K),  whereas 
IIR-SOFC  fuelled  by  methanol  presents  the  smoothest  temperature 
profile  (reformer  temperature  decreases  from  1173  to  1163  K).  It  is 
noted  that  although  the  rate  of  the  methanol  conversion  is  faster 
than  that  of  ethanol  (Fig.  3(b)  and  (c)),  the  heat  of  the  methanol 
steam  reforming  reaction  is  considerably  lower,  thus  the  temper¬ 
ature  gradient  for  IIR-SOFC  fuelled  by  methanol  is  smoother.  The 
electrical  efficiencies  achieved  from  the  IIR-SOFC  with  different 
fuels  are  compared  in  Table  3.  It  can  be  seen  that  methanol  pro¬ 
vides  the  greatest  electrical  efficiency;  nevertheless,  the  system  also 
requires  the  largest  amount  of  inlet  fuel  compared  with  other  fuels. 


Table  3 

Comparison  of  required  volumetric  flow  rate  and  electrical  efficiency  achieved  from 
IIR-SOFC  fuelled  different  fuels. 


Fuel 

Properties 

Volumetric  flow  rate 

Electrical  efficiency  (%) 

Methane 

0.15 

61.1 

Biogas 

0.18 

46.7 

Methanol 

0.22 

85.6 

Ethanol 

0.096 

32.0 
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Fig.  6.  Effect  of  operating  pressure  on  temperature  profile  of  IIR-SOFC  fuelled  by  (a)  methane,  (b)  biogas,  (c)  methanol  and  (d)  ethanol. 


Fig.  7.  Effect  of  S:C  ratio  on  temperature  profile  of  IIR-SOFC  fuelled  by  (a)  methane,  (b)  biogas,  (c)  methanol  and  (d)  ethanol. 
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Fig.  8.  Temperature  profiles  of  reformer,  fuel  and  air  channels  (Tr,  Tf  and  Ta)  for  IIR-SOFC  with  co-flow  pattern  fuelled  by  (a)  methane,  (b)  biogas,  (c)  methanol  and  (d)  ethanol. 
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Fig.  9.  Effect  of  CO2  content  in  biogas  ((a)  30%,  (b)  20%  and  (c)  10%)  on  the  temperature  profiles  of  reformer,  fuel  and  air  channels  (Tr,  Tf  and  Ta)  of  IIR-SOFC  system. 
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From  the  literature,  the  temperature  gradient  along  the  system  can 
affect  the  cell  performance  and  the  properties  of  the  ceramic  com¬ 
ponent;  the  temperature  difference  along  the  SOFC  system  should 
be  less  than  10  Kcm-1  [20].  According  to  the  present  analysis,  only 
an  IIR-SOFC  fuelled  with  methanol  can  meet  this  criterion.  Thus,  as 
another  aim  of  this  study,  several  operating  conditions,  i.e.,  operat¬ 
ing  pressure,  inlet  S:C  ratio  and  gas  flow  pattern  are  varied  in  order 
to  minimize  the  temperature  gradient  as  well  as  maximize  the  out¬ 
puts  (i.e.,  electrical  efficiency)  from  an  IIR-SOFC  fuelled  by  these 
fuels. 

3.3.  Effect  of  pressure 

According  to  the  kinetic  expressions  for  steam  reforming  of 
selected  fuels,  the  operating  pressure  plays  an  important  role  on 
gas  conversion  behaviour  as  well  as  product  selectivity.  Thus,  the 
effect  of  operating  pressure  on  the  system  performance  was  deter¬ 
mined  by  varying  the  operating  pressure  from  1  to  10  bar.  As  seen  in 
Fig.  6(a)-(d),  an  increase  in  pressure  significantly  reduces  the  tem¬ 
perature  drop  at  the  entrance  of  reformer  channel.  This  could  be 
due  to  the  effect  of  the  partial  pressure  of  gas  species  in  the  steam 
reforming  and  electrochemical  reactions  since  the  extent  of  the 
steam  reforming  rate  reaction  decreases  with  increasing  pressure 
due  to  the  thermodynamic  restrictions  and  the  kinetic  expressions 
in  use.  It  is  also  found  that  an  increase  of  operating  pressure  results 
in  a  higher  open-circuit  potential  and  consequently  lower  polar¬ 
ization,  which  leads  to  higher  electrical  efficiency.  According  to 
simulation,  the  electrical  efficiency  of  IIR-SOFC  fuelled  by  methane, 
biogas,  methanol  and  ethanol  increases  to  58.3,  44.8,  69.46,  and 
34.2%,  respectively,  when  the  operating  pressure  is  increased  from 
1  to  lObar.  Thus,  an  increase  in  operating  pressure  is  an  effective 
way  for  reducing  the  temperature  gradient,  as  well  as  enhancing 
cell  performance. 

3.4.  Effect  of  steam  to  carbon  (S/C)  ratio 

As  the  thermodynamic  properties  of  inlet  fluids  are  important 
factors  that  could  influence  the  thermal  behaviour  of  SOFC  system, 
the  effect  of  inlet  steam  content  on  the  temperature  gradient  along 
the  system  and  the  system  performance  was  investigated  by  vary¬ 
ing  the  inlet  steam  to  carbon  (S:C)  ratio  from  2.0  to  4.0.  Fig.  7(a)-(d) 
presents  the  effect  of  inlet  S:C  ratio  on  the  temperature  gradient 
along  the  reformer  channel  of  an  IIR-SOFC  fuelled  by  different  pri¬ 
mary  fuels.  It  can  be  seen  that  the  cooling  spot  reduces  when  the 
ratio  is  increased  for  all  types  of  fuel,  particularly  in  the  case  of 
ethanol.  This  improvement  could  be  caused  by  an  increase  in  heat 
accumulated  from  excess  steam  fed  into  the  system,  which  could  be 
efficiently  used  for  the  endothermic  steam  reforming  reaction.  Nev¬ 
ertheless,  it  is  also  found  that  the  electrical  efficiency  considerably 
decreases  with  increasing  inlet  S:C  ratio  (decreases  to  33.3,  29.3, 
54.6  and  15.4%  for  IIR-SOFC  fuelled  by  methane,  biogas,  methanol 
and  ethanol,  respectively,  when  the  inlet  S:C  ratio  increases  from 
2.0  to  4.0).  The  decrease  could  be  due  to  the  higher  cell  polarization 
with  increasing  inlet  S:C  ratio  and  the  higher  energy  requirement 
to  vaporize  steam  when  excess  steam  is  introduced  to  the  system 
[20]. 

3.5.  Effect  of  flow  direction 

For  typical  autothermal  application,  e.g.,  a  heat  exchanging  sys¬ 
tem,  the  flow  direction  of  exchanged  fluids  strongly  affects  the  heat 
transfer  and  reaction  behaviour  in  the  fluid  stream;  thus  the  effect 
of  the  fuel  and  oxidant  flow  direction  on  the  IIR-SOFC  performance 
was  also  considered  in  the  present  work.  In  previous  sections,  air 
flow  is  counter-flow  to  fuel  flow  in  the  fuel  channel  of  SOFC  (a  so- 
called  ‘counter-flow’  pattern).  Alternatively,  fuel  and  air  streams 


can  be  passed  in  the  same  direction,  a  so-called  “co-flow”  pattern. 
The  system  behaviour  of  with  co-flow  pattern  was  analyzed  by 
changing  mass  and  energy  balances  in  the  air  channel  along  with 
their  corresponding  boundary  conditions,  while  keeping  all  other 
operating  conditions  identical  to  those  of  the  counter-flow  pattern. 
Fig.  8(a)-(d)  shows  the  temperature  profiles  along  all  channels  of 
an  IIR-SOFC  with  co-flow  pattern.  Compared  with  the  data  in  Fig.  5, 
it  can  be  seen  that  the  flow  direction  of  fuel  gases  and  air  strongly 
affects  the  temperature  gradient  along  the  system.  The  IIR-SOFC 
with  a  co-flow  pattern  provides  smoother  temperature  distribution 
with  a  higher  average  cell  temperature,  which  leads  to  a  reduction 
in  cell  overpotentials  and  consequently  yields  a  higher  electrical 
efficiency  (increases  to  55.8, 46.6,  67.1  and  39.8%  for  methane,  bio¬ 
gas,  methanol  and  ethanol,  respectively).  The  high  efficiency  of 
an  IIR-SOFC  with  a  co-flow  pattern  is  due  to  the  good  matching 
between  the  heat  exothermically  supplied  from  electrochemical 
reaction  and  the  heat  required  for  endothermic  steam  reforming 
along  the  SOFC  system,  thus  it  is  concluded  that  an  IIR-SOFC  with 
co-flow  pattern  is  more  satisfactory  than  that  with  counter-flow 
pattern. 

3.6.  Effect  of  C02  content  in  biogas 

According  to  the  above  study,  the  electrical  efficiency  of  IIR- 
SOFC  fuelled  by  biogas  is  relatively  low  due  to  the  high  content 
of  C02  in  biogas.  It  is  well  known  that  the  presence  of  a  high  C02 
content  is  a  major  problem  with  biogas  in  terms  of  several  energy 
aspects,  e.g.,  combust  for  power  generation  or  compress  and  use 
as  fuel  in  vehicle;  thus  the  removal  of  C02  from  biogas,  which 
can  be  done  by  means  of  several  techniques,  e.g.,  water  scrubbing, 
pressure  swing  adsorption  or  membrane  technology,  has  been  pro¬ 
posed  and  widely  applied  recently.  In  the  present  work,  the  effect 
of  C02  content  in  biogas  on  system  performance  was  also  stud¬ 
ied  by  varying  the  C02  content  from  40  to  30,  20  and  10%.  It  is 
found  that  the  electrical  efficiency  increases  with  decreasing  C02 
content  (from  34.5  to  35.9, 38.2  and  39.2%,  respectively);  neverthe¬ 
less,  it  can  be  seen  from  Fig.  9(a)-(c)  that  the  content  of  C02  also 
affects  the  temperature  gradient  along  the  system.  The  tempera¬ 
ture  difference  increases  from  28.1  to  30.4,  42.8  and  47.6  K  cm-1 
when  the  C02  content  is  decreased  from  40  to  30,  20  and  10%, 
respectively,  thus  it  could  affect  the  ceramic  component  of  the 
SOFC  and  result  in  mechanical  failure  due  to  thermally  induced 
stresses. 

4.  Conclusions 

The  simulation  indicates  that,  among  the  four  types  of  fuel 
studied,  an  IIR-SOFC  fuelled  by  methanol  provides  the  greatest  per¬ 
formance  in  terms  of  electrical  efficiency  and  temperature  gradient 
along  the  system.  It  is  also  found  that  the  content  of  C02  in  biogas 
strongly  affects  the  system  performance;  by  removal  of  some  C02 
from  biogas,  the  electrical  efficiency  can  be  improved  but  a  larger 
temperature  gradient  is  expected.  By  increasing  the  operating  pres¬ 
sure,  the  system  efficiency  increases  and  the  temperature  gradient 
can  be  minimized.  The  use  of  a  high  inlet  S:C  ratio  can  reduce  the 
cooling  spot  at  the  entrance  of  the  internal  reformer  channel  but 
the  electrical  efficiency  considerably  decreases  due  to  the  higher 
cell  polarization  and  higher  energy  requirement  to  vaporize  steam 
when  excess  steam  is  introduced  to  the  system.  An  IIR-SOFC  with  a 
counter-flow  pattern  is  also  compared  with  an  IIR-SOFC  that  has  a 
co-flow  pattern.  It  is  found  that  a  co-flow  pattern  provides  higher 
voltage  a  and  smoother  temperature  gradient  along  the  system  due 
to  better  matching  between  the  heat  supplied  from  the  electro¬ 
chemical  reaction  and  the  heat  required  for  the  steam  reforming 
reaction. 
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